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Quiescent X-ray Emission from an evolved Brown Dwarf ? 

B. Stelzeri'2 
stelzerSastropa . unipa . it 

ABSTRACT 

I report on the X-ray detection of Gl 569 Bab. During a 25ksec Chandra 
observation the binary brown dwarf is for the first time spatially separated in X- 
rays from the flare star primary Gl 569 A. Companionship to Gl 569 A constrains 
the age of the brown dwarf pair to ~ 300 — 800 Myr. The observation presented 
here is only the second X-ray detection of an evolved brown dwarf. About half of 
the observing time is dominated by a large flare on Gl 569 Bab, the remainder is 
characterized by weak and non-variable emission just above the detection limit. 
This emission - if not related to the afterglow of the flare - represents the first 
detection of a quiescent corona on a brown dwarf, representing an important 
piece in the puzzle of dynamos in the sub-stellar regime. 

Subject headings: X-rays: stars. Stars: coronae, flare, brown-dwarfs, individual 
(G1569Bab) 

1. Introduction 



Late-type stars have long been known to display signatures of magnetic activity evi- 
dencing solar-like dynamo action (Rosner, Golub, & Vaiana 1985). The picture is less clear 
for the coolest dwarf stars at the bottom of the main-sequence (MS) and in the substellar 
regime. These objects are fully convective and can not drive a solar-like afi-dynamo. The 
efficiency of alternative field generating mechanisms, e.g. an a^-dynamo (Radler et al. 1990) 
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or a turbulent dynamo (Durney, De Young, & Roxburgh 1993), are both difficult to pre- 
dict and observationally poorly constrained. Systematic investigations of chromospheric Ha 
activity in very low-mass (VLM) stars and brown dwarfs (BDs) in the solar neighborhood 
(Gizis et al. 2000; Mohanty & Basri 2003) have shown that Ha emission reaches a maxi- 
mum at spectral type M7, and fades off for the latest M dwarfs suggesting a change in the 
driving dynamo mechanism. X-ray emission is a complementary activity indicator, probing 
the effects of the magnetic processes in the outermost and hottest part of the atmosphere, 
the corona. Beyond spectral type M6 the X-ray regime is widely unexplored. Only few late 
M dwarfs have been observed with Chandra and/or XMM-Newton: APMPM2354 (Scholz et 
al. 2004), vBlO (Fleming, Giampapa, & Garza 2003), LP944-20 (Rutledge et al. 2000). 
With the exception of vB 10 all X-ray detections of M dwarfs cooler than spectral type M7 
are ascribed to flares, i.e. outbursts on a time-scale of hours that arise from reconnection of 
magnetic field lines. It is unclear to date whether ultra-cool dwarfs can sustain persistent 
X-ray emission, as is typically observed on higher-mass stars. In particular, a 'quiescent' 
corona could not be established for any BD so far. 

While the observational material on X-ray activity of evolved ultra-cool dwarf stars 
and BDs is scarce. X-ray emission has lately been revealed from a substantial number of 
~ 1 — 5 Myr-old VLM pre-MS stars and BDs in various star forming regions (Neuhauser 
& Comeron 1998; Preibisch & Zinnecker 2002; Mokler & Stelzer 2002; Feigelson et al. 
2002; Stelzer, Micela, & Neuhauser 2004), for one BD in the ~ 10 - 30 Myr-old TWHya 
association (Tsuboi et al. 2003), one member of Upper Sco (Bouy 2004), and one BD in the 
Pleiades (Briggs & Pye 2004). BDs are not able to fuse hydrogen and subsequently must 
cool down and become fainter as they age. Therefore, their detection at young ages but 
non-detection in more evolved stages could be due to: (A) a drop of the X-ray luminosity 
as a result of the decreased bolomctric luminosity. A decline of with age is expected 
if the Lx/Lboi correlation typical for the higher- mass late-type stars holds also in the sub- 
stellar regime; (B) an effect of the atmospheric temperature. Based on the absence or low 
level of Ha emission seen in most field L dwarfs Mohanty et al. (2002) have argued that 
the dynamo may shut off below a critical temperature because the atmosphere becomes too 
neutral to provide substantial coupling between matter and magnetic field. But due to a lack 
of systematic observations this hypothesis has not been tested so far in the X-ray regime, 
nor has a firm relation between X-ray and Ha emission been established. 

For a study of the evolution of X-ray activity on sub-stellar objects as they evolve 
away from the youngest ages BD companions to MS stars are favorable targets. The age 
determination in the coolest part of the HR diagram is subject to considerable uncertainty 
because of a degeneracy in the mass-age plot and a uncertainty of theoretical evolutionary 
models for VLM stars and BDs. But ages of MS stars can be estimated by a number of 
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methods (Lachaume et al. 1999), and companions are expected to be coeval. 

A low-mass companion to the nearby (9.8 pc; Dahn et al. 2002) dM2e star G1569 
was identified and confirmed as a common proper motion pair by Forrest, Skrutskie, & 
Shure (1988). Martm et al. (2000) found that the companion itself is a binary of 0.1" 
separation. Various methods have been used to constrain the age of G1569A: Ha and 
X-ray activity (Henry & Kirkpatrick 1990; Pallavicini, Tagliaferri, & Stella 1990), the 
non-dctcction of Lithium absorption (Magazzii, Martin, & Rcbolo 1993), and its space 
velocity suggesting membership to the UMa association or a new super-cluster (Chereul, 
Crczc, & Bienaymc 1999). All these methods place it at an age between 300 — 800 Myr, 
consistent with the position of the Gl 569 Bab binary in the HR diagram (Lane et al. 2001). 
The dynamical masses of G1569Ba and G1569Bb were derived by Zapatero Osorio et al. 
(2004) from a combination of radial velocity measurements and astrometry, yielding 99% 
confidence intervals of Mga = 0.055-0.087 M© and Mgh = 0.034-0.070 M©. The lower-mass 
component is thus the first object confirmed as substellar independent of models. In this 
Letter I present the Chandra detection of Gl 569 Bab, including a large flare and probably 
the first detection of quiescent X-ray emission from an evolved BD. 



2. Observations 

The G1569 system was observed with the Advanced CCD Imaging Spectrometer (ACIS), 
using the S3 chip in imaging mode for a total exposure time of 25 ks. The data analysis 
was carried out using the CIAO software package^ version 2.3 in combination with the 
calibration database (CALDB) version 2.21. I started the analysis with the level 1 events 
file provided by the pipeline processing at the Chandra X-ray Center (CXC). In the process 
of converting the level 1 events file to a level 2 events file for each of the observations the 
following steps were performed: Charge transfer inefficiency (CTI) correction, removal of the 
pixel randomization, filtering for event grades (retaining the standard ASCA grades 0, 2, 3, 
4, and 6), and apphcation of the standard good time interval (GTI) file. The events file was 
also checked for any systematic aspect offsets using CIAO software, but none were present. 



^CIAO is made available by the CXC and can be downloaded from 
http: / / cxc.liarvard.edu/ ciao/download-ciao-reg.html 
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3. Data Analysis and Results 

The expected position of G1569 A was found by translating its Hipparcos position to the 
time of the Chandra observation with help of the proper motion given by King et al. (2003). 
The position of Gl 568 Bab relative to the primary was obtained using the separation and 
P. A. of Lane et al. (2001). Source detection was carried out with the wavdetect algorithm 
on an image of 50 x 50 pixels length (1 pixel = 0.492") centered on the computed position 
of the primary using wavelet scales between 1 and 8 in steps of -\/2. Two X-ray sources are 
found which are unambiguously identified with G1569A (offset between optical and X-ray 
position Aox = 0.28") and Gl 569 Bab (Aqx = 0.45"), respectively; see Fig. 5. Using the 
position of the X-ray source corresponding to G1569A as the reference, the second X-ray 
source is characterized by an angular separation of (5.02 ± 0.06)" and a position angle of 
(31.7 ± 0.6)°. These values are in reasonable agreement with those published by Lane et al. 
(2001) based on near-IR images, the difference possibly being due to orbital motion of the 
BD pair around the primary in the time interval of ~ 3.4yr. 

I extracted the source photons of Gl 569 Bab from a circular area of 1.5" radius centered 
on the position of the corresponding X-ray source. According to a simulation of the point 
spread function (PSF) with mkpsf this area contains about 95 % of the source photons assum- 
ing a photon energy of 1.5 keV. To eliminate contamination by counts from the wings of the 
PSF of the X-ray bright primary the background was computed within an annulus centered 
on Gl 569 A and at the distance of the BD binary, but excluding the source extraction area 
of Gl 569 Bab. After correcting for the missing 5% of photons the total number of source 
counts is 250. 

The primary G1569 A suffers from pile-up, such that the flux is underestimated and the 
shape of the lightcurve and the spectrum are distorted when source photons are extracted 
from a circular area centered on the X-ray source. At the expense of reduced statistics the 
negative effects of pile-up can be avoided if photons from the core of the PSF are ignored. To 
determine the optimum inner radius for an annular source extraction region I extracted 
events from a series of annuli with rin = 0.5".... 2. 25". The outer radius was fixed at 3.5" 
to exclude photons from the BD pair. With help of the ChaRT^ and MARX^ simulators 
an individual ar/ was constructed for each of the annuli to correct the effective area for the 



^Information on the usage of the Chandra Ray Tracer (ChaRT) can be found at 
http: / / cxc.harvard.edu/ chart/ 

^The Model of AXAF Response to X-rays (MARX) is provided by the MIT/CXC (see 
http://space.niit.edu/CXC/MARX/) 
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missing part of the PSF^. Then I examined how the observed luminosity, derived from the 
photon spectrum, changes depending on Tin- The optimized extraction area was then defined 
as the annulus where Lx reaches a plateau, rin = 1.25". 

The X-ray hghtcurves of Gl 569 Bab and G1569A are shown in Fig. 5. Both the BD 
pair and the M2 star underwent a spectacular flare. Note that the relative amplitude of the 
two lightcurves does not reflect the actual strength of the outbursts because - as explained 
above - for G1569 A photons from only a small fraction of the PSF (4.5 % for photon energy 
of 1.4967 keV) could be considered in the source extraction region. The strength of pileup 
depends on the changing incident photon flux. Therefore the exact time evolution of the 
intensity of G1569 A is difficult to reconstruct. The X-ray emission of Gl 569 A is not the main 
scope of this Letter, and will be discussed elsewhere. In the following I focus on Gl 569 Bab. 



4. X-ray Properties of Gl 569 Bab 

The beginning of the observation is dominated by a strong flare on the BD binary, 
with --^ 93% of all source photons being concentrated in the flrst 12ksec of the observation. 
In the remainder of the observing time a total of 23 photons were collected in the source 
area, of which about 6 are expected to be background events, mostly contaminants from 
G1569 A. This corresponds to a 'quiescent' source count rate of (1.2±0.3) 10~^cps. Assuming 
a 1-T plasma without absorption yields logLg = 25.8erg/s. This value is insensitive to 
temperature within a typical range of 0.3...1.2 kcV. The peak count rate diiring the flare 
is ~ 0.0625 cps, amounting to a peak luminosity of logLp = 27.5erg/s using the same 
assumptions for the spectral shape as for the estimate of the 'post-flare' luminosity. A more 
reliable estimate for the X-ray luminosity can be obtained from an analysis of the X-ray 
spectrum. It provides information on the time-averaged luminosity, which is dominated by 
the flare. A statistically acceptable interpretation in terms of thermal emission requires 
a model with a minimum of two temperatures {kT^ = O.SStoH keV, kT2 = 0.94l°:| keV, 
EM1/EM2 — 1.2lo;g). Discrepancies between model and data suggest that absorption may 
not be neghgible (see Fig. 5). However, if an absorption term is added, the fit results 
in very high emission measure for the soft component. A correlation in the parameter 
space between A^h and EMi is typical for low-sensitivity, low-resolution spectra. The model 
without absorption seems to be more realistic, because (i) no visual extinction has been 
reported for G1569, and (ii) a dominant soft component is unexpected in flares. 



see http://www.astro.psu.edu/users/tsujiniot/arfcorr.htmlfor more information on the procedure 
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Table 1: X-ray luminosity and L^/ L^^i ratio for Gl 569 Bab. 





Quies 


Peak 


Average 


log -Lx 


[erg/s] 25.8 


27.5 


26.8 


log (-Lx/i-boOBa 


-4.4 


-2.7 


-3.4 


log (-Lx/i-boOBb 


-4.2 


-2.5 


-3.2 
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A widespread used measure for coronal activity is the fraction of X-ray luminosity rela- 
tive to the bolometric luminosity. The two components in the BD binary have log Lea [erg/s] = 
30.2 and logLeb [erg/s] = 30.0 (Gorlova et al. 2003). Table 1 presents a summary of X- 
ray luminosities and Lx/Lboi ratios for Gl 569 Bab derived from this Chandra observation. 
Lx/Lboi was evaluated attributing all observed X-rays to either of the unresolved compo- 
nents. Recall that, 'Quies' and 'Peak' values are estimated from count rates using model 
assumptions. 'Average' denotes the luminosity derived from the analysis of the time averaged 
spectrum. 



5. Discussion 

I report on a huge X-ray flare and the possible detection of quiescent emission on 
G1569 Bab. The present 3a limits on the individual masses do not exclude that the higher- 
mass component is a star at the borderline to the substellar regime, but most evidence points 
at a BD-BD binary (see discussion in Zapatero Osorio et al. 2004). Independent of whether 
the X-ray source is substellar or not, the emitter is one of only two ultra-cool dwarfs older 
than ~ 100 Myr for which X-ray emission could be detected so far. Therefore, Gl 569 Bab 
constitutes an important link between activity in star forming regions, associations, and 
open clusters, and the more evolved field stars and BDs. Gl 569 Bab are spectroscopically 
confirmed companions to a MS star. This eliminates two important parameters (the effec- 
tive temperature and the age) from the list of unknowns, whose influence on the efficiency 
of (sub-)stellar dynamos has remained elusive so far. 

Fig. 5 shows the Lx/Lboi vs. spectral type diagram for VLM field dwarfs. Gl 569 Bab 
adds to a handful of late-M type objects that have shown fiares, with emission levels near or 
above the canonical 'saturation level' of 10~^. In the second half of the Chandra observation, 
after the flare decayed, Gl 569 Bab is detected as a weak X-ray source. If the interpretation 
of this emission as representative for the quiescent corona is correct it constitutes the flrst 
detection of persistent X-ray activity on a presumed BD. The sharp decline in the quiescent 
coronal activity for ultra-cool dwarfs conjectured by Fleming, Giampapa, & Garza (2003) 
must then be questioned. As seen from the upper limits in Fig. 5, indeed, most previous 
observations of late M dwarfs were not sensitive enough to sample the range expected from 
an extrapolation of the Lx/Lboi values of the M0...M6 dwarfs. However, the non-detection 
of any photon from Gl 569 Bab in the ~ 2.2 ks before the flare casts some doubt on this 
view. It is observed occasionally that stars do not go back into their pre-flare quiescent state 
immediately after an outburst, thus we may see some kind of 'afterglow' of the flare rather 
than true quiescent emission. On the other hand, based on the post-flare emission level a 



-8- 



total of 2.6 source photons are expected prior to the flare, consistent with the observation 
of zero photons at the 2(7 hmit. 

The detection of a strong X-ray flare on a comparatively evolved BD demonstrates that 
coronae of VLM objects can remain powerful beyond the youngest ages. Fig. 5 puts the 
Chandra observation of Gl 569 Bab in context with observations of younger BDs, and the 
only previous X-ray detection of a field BD, LP 944-20. A factor of 4 — 5 more photons 
have been collected for Gl 569 Bab with respect to all other late M dwarfs, allowing for the 
first time a meaningful spectral analysis of such an object. The spectrum shows significant 
emission above IkcV, but a similarly strong cool component of ~ 3.2 MK. Low coronal 
temperatures seem to be characteristic for all but the youngest BDs. This is all the more 
remarkable as temperatures are known to rise in stellar flares. The flare temperatures of 
the evolved BDs (Gl 569 Bab and LP 944-20) are similar to the quiescent temperature of the 
middle-aged BD TWA-5B on the one hand and to the quiescent temperature of the evolved 
star vB 10 on the other hand. Upcoming observations with Chandra sampling the age and 
temperature space shall constrain the influence of these parameters on coronal emission. 

I would like to thank G. Micela and E. Flaccomio for fruitful discussions and careful 
reading of the manuscript. 
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Fig. 1. — Chandra ACIS-S image of the G1569 system. The size of the image is 50 x 
50 pixels. X-points denote the optical/IR position of G1569A and Gl 569 Bab. EUipses 
mark the corresponding X-ray sources. 
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Fig. 2. — Source and background lightcurves of Gl 569 Bab and G1569A. 
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Fig. 3. — X-ray spectrum of Gl 569 Bab, 2-T APEC model, and residuals. 
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Fig. 4. — Ratio of L^/ Lboi versus spectral type for VLM field dwarfs, filled circles - M dwarfs 
detected during the RASS (see Hiinsch et al. 1999); asterisks - flares, squares - quiescent 
emission or upper limits (data from F93, F95, N99, FOO, ROO, SL02, F03, H04, S04); large 
asterisk and square - new data for Gl 569 Bab. For clarity VLM dwarfs with the same spectral 
type have a small offset with respect to each other on the horizontal axis. 
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Fig. 5. — Ratio of Lx/Lboi versus age for BDs: triangles, diamonds, open circles - bona-fide 
BDs in star forming regions (spectral types > M7) with a random one order of magnitude 
spread in age centered on 0.5 Myr for p Oph (X-ray data from 101), 1.5 Myr for IC 348 (PZ02), 
and 2 Myr for Chal (St04); flaring BDs are displayed as an asterisk, and those with accretion 
signatures in the H« profile as x-points. References for evolved BDs: DENIS J1556 (B04), 
TWA-5B (T03), 2MW-J1207 (GB04), Roquc 14 and four undetected BDs in the Pleiades 
(BP04), LP944-20 (ROO), and Gl 569 Bab (this paper). 



